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We use QCD sum rules, in the framework of the Heavy Quark Effective The- 
ory, to calculate the universal form factor ^^(y) parameterizing the semilep- 
t^-. | tonic transitions B — ► DqIv, B — » D*£i>, where Do and D\ are the members 

of the excited charmed doublet with J p = (0]y 2 , l]y 2 )- We include two- loop 
^-j- ■ corrections in the perturbative contribution to the sum rule, and present a 

■ complete next-to- leading order result. As a preliminary part of our analysis 

00 ■ 

we also compute, up to order a s , the leptonic constant F + of the doublet Dq, 
D*. Finally, we discuss the phenomenological implications of this calculation. 
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I. INTRODUCTION 



The application of the heavy quark flavor and spin symmetry, valid in QCD in the 
infinite heavy quark mass limit together with the heavy quark effective field theory 

(HQET) 0-0, has led to a dramatic progress towards a model-independent description 
of the spectroscopy and the decays of hadrons containing a single heavy quark Q (Q = 
c,b). An outstanding result of the theory concerns the description of the exclusive B — > 
Div and B — > D*iv semileptonic decays, in the limit vhq — > oo, in terms of just one 
nonperturbative, universal form factor (the Isgur-Wise function £), normalized to unity 
at maximum momentum transfer to the lepton pair. Other distinctive examples are the 
relations between the beauty meson leptonic constants and the beauty meson semileptonic 
transition amplitudes to light mesons at zero recoil, with the analogous charmed meson ones, 
obtained employing general dimensional scaling rules. 

Corrections of order KqcD/ m Q (and higher powers) to the leading term can be system- 
atically analyzed in HQET in terms of a reduced number of hadronic, universal parameters, 
with a remarkable simplification of the analysis. However, in the applications of HQET the 
effects of non-perturbative strong interactions can be estimated only in the framework of 
some non-perturbative theoretical approach. In this regard, particularly fruitful has been 
the application of sum rules || formulated in the framework of HQET [^]. This method 
is genuinely field theoretical and based on first principles, and relates the hadronic ob- 
servables to QCD fundamental parameters via the Operator Product Expansion (OPE) of 
suitable Green's functions. Such an expansion involves perturbative contributions as well 
as non-perturbative quark and gluon vacuum condensates. In particular, a s corrections to 
the coefficients of the OPE can be computed order by order in perturbation theory, and 
therefore they can be systematically taken into account. 

A critical aspect of the sum rule calculations in HQET is represented by the size of 
non- leading terms, such as the 1 / tuq corrections and the a s corrections in the perturbative 
expansion of the OPE. For example, the predictions for the leptonic constants of qQ pseu- 
doscalar mesons are affected by considerably large next-to- leading corrections in « s |8|-10 



also AqcD/ m Q corrections are non- negligible in the case of the D meson, an effect confirmed 
by lattice QCD analyses [[□] . 



Conversely, in the HQET QCD sum rule calculation of the Isgur-Wise function, the next- 



to-leading order a s corrections turn out to be small and well under control [12,13], and the 



same is true for Aqcd/it^q corrections [ 14| , specially near the zero recoil point where the 



normalization of the universal form factor is protected by the heavy quark symmetry. This 
has allowed a drastic reduction of the theoretical uncertainty in the determination of the 
CKM matrix element V c b ITS . 



It is worth analyzing other cases analogous to the determination of the Isgur-Wise form 
factor £, and we present here a HQET sum rule calculation of the universal form factor 
governing the semileptonic B meson decays into the + , 1 + charmed excited states, up to 
next-to- leading order in a s and to leading order in the heavy quark expansion mq — > oo. 
These higher-lying charmed states correspond to the L = 1 orbital excitations in the non- 
relativistic constituent quark model. Besides their theoretical relevance to HQET [16|], in 
particular to the aspects of the QCD sum rule calculation mentioned above, such B — > D** 
semileptonic transitions (D** is the generic L = 1 charmed state) have numerous additional 
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points of physical interest. 

Indeed, in principle these decay modes may account for a sizeable fraction of semileptonic 
.B-decays, and consequently they represent a well-defined set of corrections to the theoretical 
prediction that, in the limit mq — > oo and under the condition (m^ — m c ) / {m^+m^j — > (the 
so-called small- velocity limit), the total semileptonic B — > X c decay rate should be saturated 
by the B —> D and B — > D* modes ||. Moreover, the shape of the inclusive differential 
decay distribution in the lepton energy could reflect contributions from the B — ► D** modes. 

Another important result, relevant both to phenomenology and to the critical tests of 
HQET, is the relation of the B — > D** form factors at zero recoil to the slope of the B — > 
Isgur-Wise function, through the Bjorken sum rule |T7|]. Of similar interest for HQET is the 
test of the upper bound on such universal form factors at zero recoil, involving the heavy 
meson "binding energy" and the D** — D mass splittings, that is the analog of the "optical" 
sum rule for dipole scattering of light in atomic physics |T8l , |l9 . 



Moreover, the KqcD/ m Q corrections can have a role for S-decay modes into excited 
charmed states, that mostly occur near the zero recoil point where the corresponding transi- 
tion matrix elements vanish. The shape of the lepton energy spectrum near such kinematical 
point including the 1 / tuq corrections, that in HQET can be predicted in terms of the Isgur- 
Wise function and mesons mass splittings, represents an important test of the theory |f2U . 



Continuing with the aspects justifying the interest for B — > D**, let us notice that 
the investigation of the semileptonic B transitions to excited charm states is an important 
preliminary study for the theoretical analysis of the production of such states in nonleptonic 
B decays as well as for the identification of additional decay modes (such as D^D^tt) 
suitable for the investigation of CP violating effects at B factories E3]. 

Finally, as a byproduct of the QCD sum rule calculation, theoretical predictions about 
the yet unobserved D** meson masses can be obtained, that are obviously interesting per 
se. 

In the following we present a complete next-to-leading order evaluation of the 5-meson 
semileptonic transition to the scalar charmed state by QCD sum rules, at the leading order 
of m^ c — > oo. In Sect. II we report the main aspects of the spectroscopy and decays of L = 1 
qQ mesons, together with the definition of the universal form factor Tuzty). The various 
steps of the QCD sum rule determination of such form factor, within HQET, are collected 
in Sect. Ill and V-VII, together with the analysis, in Sect. IV, of the leptonic constant F + 
of the doublet Dq,D\. In Sect. VIII the phenomenological implications of our calculation 
are presented, together with the conclusions. 



II. POSITIVE PARITY HEAVY-LIGHT MESONS 

In the infinite heavy quark mass limit the spectrosopy of hadrons containing one heavy 
quark Q is greatly simplified, due to the decoupling of the heavy quark spin sq from the 
angular momentum of the light degrees of freedom (quarks and gluons) sg = J — sq. This 
allows a classification of such hadronic states by J and by S£, so that hadrons corresponding 
to the same sg belong to degenerate doublets. In the case of qQ mesons, the low-lying states 
with sf = g correspond to the pseudoscalar 0~ and vector 1~ mesons (B, B*; D, D*), the s- 
wave states of the constituent quark model. The four states corresponding to orbital angular 
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momentum L 



(^1/2 > I1/2) anc ^ J P = (-*-3/2> 2 3 / 2 )> 
which differ by the values sf = | + and sf = | + , respectively. The states 0t 2 and 2^ 
correspond to the scalar and spin-two mesons of the quark model; the relation between 
HQET and quark model L 



1 can be classified in two doublets: J p 
and "'' 



1, 1 + states can easily be derived [23 



The charmed 2^ 2 state has been experimentally observed and denoted as the 1)2(2460) 
meson, with m D * = 2458.9 ± 2.0 MeV, T D * = 23 ± 5 MeV and m D * = 2459 ± 4 MeV, 
Yd* = 25^7 MeV for the neutral and charged states, respectively |[24| . The HQET state 
l+ /2 can be identified with L>i(2420), with m Dl = 2422.2 ± 1.8 MeV and T Dl = 18.9 

MeV " 



+4.6 
3.5 



24], even though a 1, 



Ll y 2 component can be contained in such physical state due to the 
mixing allowed for the finite value of the charm quark mass. [] 



Both the states 2^ 2 and 



lg/ 2 decay to hadrons by d— wave transitions, which explains 



their narrow width; the strong coupling constant governing their two-body decays can be 
determined using experimental information [26]. On the other hand, the sf = | + doublet 
(D , Z?*) has not been observed yet. The strong decays of such states occur through s- 
wave transitions, with expected larger widths than in the case of the doublet | + . Indeed, 
analyses of the coupling constant governing the two-body hadronic transitions by QCD sum 
rules predict T(D° -> D + tj-) ~ 180 MeV and T(Df -> D* + tt~) ~ 165 MeV [|7| . Estimates 
of the mixing angle a between D* and Di give a ~ 16° [2"7,2"8fl. 

The matrix elements of the semileptonic B — > DqCv and B — > D\lv transitions can be 
parameterized in terms of six form factors: 



< D (v')\crr^ 5 b\B(v) 



> 



y/m B m Do 
<D\{v', e)\crf^l- lh )b\B{v) > 



9+{v + v) fl + g- 



9vi^ + e* ■ v [g V2 v^ + gvsV^] - % ^W^^V 7 , (2.1) 



where v and v' are four-velocities and e is the D\ polarization vector. The form factors gi 
depend on the variable y = v ■ v', which is directly related to the momentum transfer to the 
lepton pair. In terms of such form factors the semileptonic differential decay rates can be 
expressed as 



dT 

dy 



dT 



(B -> D £u) 



dy 



G 2 F \V cb \ 2 m% 
48tt 3 

G 2 F \V cb \ 2 m % 
48tt 3 



V 7Tl B ' LV 



v^T(^) 3 { [(» - + (:</ - D(jv, + 



1 + 



m B 



9+-(l 



m B 1 J 



m B 



m B 



m B 



+ 2[l-2^+(^)>U(!/ 2 -l)sy} • 
L m B v m B / J 1 A > 



■ m D *^2 



(2.2) 



The heavy quark spin symmetry allows to relate the form factors gi{y) in ( [2.1|) to a 
single function Ti/ 2 (y) [jTB] through short-distance coefficients, perturbatively calculable, 



1 There is also experimental evidence of beauty s\ 



3 + 
2 



states [25|. 
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which depend on the heavy quark masses m&, m c , on y and on a mass-scale fi, and connect 
the QCD vector and axial vector currents to the HQET currents. At the next-to-leading 
logarithmic approximation in a s and in the infinite heavy quark mass limit, the relations 
between and 71/2 are given by 

9+(y) + 9-(v) = "2 (CUV) + (y- l)Cf (//)) r 1/2 (y, /i) 

9+(y) - 9-(v) = 2 (c?(/0 - (y - l)Cf 0*)) r 1/2 (y, 
^Vi(y) = 2(y- 1) Ci(/x) T 1/2 (y,fi) 
9v 2 (y) = -2 C 2 (/i) T 1/2 (y,n) 
9v z {y) = -2 (Ci(/i) + C 3 (/i)J T 1/2 (y,n) 

g A (y) = -2C!(n)T l/2 (y,fi) . (2.3) 

The //-dependence is the same for all the functions CI , where Ci refer to the vector current 
and Cf to the axial one; therefore, one can extract such a dependence by writing fTj 

Cf V) = % s \m b ,m B ,y) K hh (y,fi) (2.4) 

where K hh = Mjx)] {l - ^Z M (y)}, with a hh = § 7 (y), 

T(2/) = ^k(y)-l] (2.5) 



and r(y) = lri ( y+ ^/ y J2 ; j s related to the velocity-dependent anomalous dimension of the 



heavy-heavy 6 — > c current in HQET [Q. The coefficient Z^, derived in |f7|j9|j29[1 has the 
expansion, for rif = 3: 

,_ /752 87r 2 \. 1N / 368 4tt 2 \. 1N2 , o e . 



The next-to-leading- log expression of the coefficient functions C^ can be found in [0. At 



the leading-log approximation they simply read: C\ = Cf = 2 (m ) [ a s(n^c)] ahh , with 
ai = — Jr, the coefficients and C^ 5 "* being zero. 



°J r 



From (|2.3|) and (|2.4|) it follows that the product Khh(y, fi^Ti/^y, a0 does not depend on 
// , and therefore a renormalization-group invariant form factor can be defined 

T^(y) = K hh (y,riT 1/2 (y,v); (2.7) 

hence, in ( |2.3|) one can substitute Cj; 5 \n) with C'f 5 '' and Ty 2 (y,ii) with T[J 2 (y). 

Analogous relations hold for the eight form factors parameterizing the matrix elements 
of B — > D\lv and i? — > D\iv\ in this case the heavy quark symmetry allows to relate 
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them to another universal function r 3 / 2 (y) fll6| . The main difference with respect to the 
Isgur-Wise form factor £(y) is that one cannot invoke symmetry arguments to predict the 
normalization of both Ti/ 2 (y) and t 3 / 2 (?/), and therefore a calculation of the form factors 
in the whole kinematical range is required. For B — > (D Q , Dl)£u the physical range for 
the variable y is restricted between y = 1 and y = 1.309 — 1.326, taking into account the 
values for the mass of D , D\ {mr> ,D* = 2.40 — 2.45 GeV). Consequently, for practical 
purposes one can adopt a linear parameterization of the form factor 71/2(2/): 7i/ 2 (2/) — 
Ti/ 2 (1)(1 — Py 2 (y ~ !))• This parameterization allows an easier comparison between the 
predictions of different approaches. 



A determination of 71/2(2/) by QCD sum rules at 0(a s = 0) was carried out in P0 
starting from finite values of the beauty and charm quark masses and the performing the 
limit m&, m c — > 00. The obtained result can be summarized by the values 7~i/ 2 (l) ~ 0.25 and 
~ 0.4, which imply quite small values of the branching ratios of B — > (D , Dl)£u. 

Other determinations of 71/2(2/) have appeared in the literature, employing various ver- 
sions of the constituent quark model |pT|-|35||. The results range in a quite large interval, 
71/2(1) = 0.06 — 0.40 and p^ 2 = 0.7 — 1.0, and critically depend on the peculiar features of 
the models employed in the numerical calculation. 

As for T 3 / 2 (y), a QCD sum rule analysis to the leading order in a s = |31J gives r 3 / 2 (l) ~ 
0.28 and p\, 2 ~ 0.9. Quark model results, on the other hand, give predictions in the range 
r 3 / 2 (l) ~ 0.31 — 0.66 and p\j 2 — 1.4 — 2.8 f3l|-|3"5| . We do not consider here the problem 



of the role of radiative corrections to the function r 3 / 2 , but limit our analysis to the case of 
T1/2, for which a number of interesting information can be worked out. 

In the next Section we briefly outline the basic points of the QCD sum rule method, as 



needed for the extension of the calculation of |30| to the next-to-leading order in a 



III. FORM FACTOR r 1/2 FROM QCD SUM RULES IN HQET 

Following 0, the determination of the universal function Ti/ 2 (y) by QCD sum rules in 
HQET is based on the analysis of the three-point correlator 

rL>,u/,y) =% 2 I dx dzeW x ~ k ^ < 0|T[j/(x), I M (0), J^}\0 > 



= i(v-v') li IL(w,u/,y) , (3.1) 

where = ^q/T^Ts^q is the b — > c weak axial current, and Jf = qhg, and J5 = qi^hg 
represent local interpolating currents of the scalar (Do) and pseudoscalar (B) mesons in 
eq. ( |2.1| ) represented in terms of HQET Kq fields and light quark q fields. The variables 



k, k' are "residual" momenta, obtained by the expansion of the heavy meson momenta in 
terms of the four- velocities: P = vtlqv + k, P' = itlqiv' + k'\ they are O(Aqqd), and remain 
finite in the heavy quark limit. 

Using the analyticity of H(w, 00', y) in the variables u = 2v ■ k and u' = 2v' ■ k! at fixed 
y, one can represent the correlator fl3.1|) by a double dispersion relation of the form 



"For an extensive discussion of qQ meson interpolating currents in HQET see, e.g., [36| 
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n(u,u',y)= \dvdv'- Pi ^') y) - , (3.2) 

J [y — uj — le) \y' — uj' — it) 

apart from possible subtraction terms. The correlator H(uj,uj',y) receives contributions 
from poles located at positive real values of uj and uj' , corresponding to the physical single 
particle hadronic states in the spectral function p[y, v' ,y). The lowest-lying contribution is 
represented by the 0~,0 + qQ and qQ' states, i.e. B and D . This contribution introduces 
the form factor T\/2 through the relation: 

U ^ W ' V) = (2A-a;-ie)(2A+-a/-i e ) ' (3 ' 3) 

where \i is a renormalization scale and F(n), F + (p) are the couplings of the pseudoscalar 
and scalar interpolating currents to the 0~ and + states, respectively, in the heavy quark 
theory: 

< 0\J%\B(v) > = F(fi) (3.4) 
<0\J v s \B (v) > =F + (ii) . (3.5) 

F(fj,) and F + (fi) are scale-dependent low energy HQET parameters, which do not depend 
on the heavy quark mass mq or rrig. In particular, F(fi) is related to the B- meson leptonic 
decay constant The mass parameters A and A + identify the position of the poles in uj 
and uj', and can be interpreted as binding energies of the 0~ and + states: A = M B — m b , 
A+ = M Do - m c . 

The higher states contributions to p(u, u',y) can be taken into account by a QCD con- 
tinuum starting at some thresholds v c and v' c , and are modeled by the asymptotic freedom, 
perturbative spectral function fp eri {y, v\y) according to the quark-hadron duality assump- 
tion. Here, p pert is the absorptive part of the perturbative quark-triangle diagrams, with 
two heavy quark lines corresponding to the weak b — > c vertex and one light quark line 
connecting the heavy meson interpolating current vertices in ( |3.1|) . At the next-to- leading 
order in a s , all possible internal gluon lines in such triangle diagrams must be considered. 

Therefore, for the dispersive representation ( p.2|) in terms of hadronic intermediate states 
one assumes the ansatz 

U(u, u', y) = Tl P oie{u, uj', y) + U continuum (uj , uj', y) (3.6) 

where, for simplicity, the dependence of the continuum contribution on the thresholds v c 
and v' c has been omitted. 

The correlator IT(u;, uj' , y) can be expressed in QCD in the Euclidean region, i.e. for large 
negative values of ui and uj' , in terms of perturbative and nonperturbative contributions: 

n(w, uj', y)= [ dvdv'- P'^^y) + U n P{uJ: J, y) . (3.7) 

J [y — ui — it) [y — uj' — it) 

In ( |3.7| ) Il np represents the series of power corrections in the "small" - and variables, 
determined by quark and gluon vacuum condensates ordered by increasing dimension. These 
"universal" QCD parameters account for general properties of the nonperturbative strong 
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interactions, for which asymptotic freedom cannot be applied. The lowest dimensional ones 
can be obtained from independent theoretical sources, or fitted from other applications 
of QCD sum rules to cases where the hadronic dispersive contribution is particularly well 
known. In practice, since the higher dimensional condensates are not known, one truncates 
the power series and a posteriori verifies the validity of such an approximation. In our 
application we shall include the dimension three quark condensate and the dimension five 
quark-gluon mixed condensate, which are known rather reliably; we neglect the contribution 
of the gluon condensate, which always turns out to be numerically small in the analysis of 
heavy-light meson systems. The calculation of fP ert and H np at the next-to-leading order a s 
is described in the sequel. 

The QCD sum rule for T\/2 is finally obtained by imposing that the two representations of 
U(u,u',y), namely the QCD representation ( |3.7| ) and the pole-plus-continuum ansatz fl3.6|), 
match in a suitable range of Euclidean values of uj and uj'. 

A double Borel transform in the variables uj and uj 1 

—j gM _ n m — ( ^_\ n ( n ^ u _^ — oo, r = fixed) (3.8) 

t (n — 1)! duj n 

(and similar for B T >) is applied to "optimize" the sum rule. As a matter of fact, this 
operation has two effects. The first one consists in factorially improving the convergence 
of the nonperturbative series, justifying the truncation procedure; the second effect is to 
enhance the role of the lowest-lying meson states while minimizing that of the model for the 
hadron continuum. The a priori undetermined mass parameters r and r' must be chosen 
in a suitable range of values, in the present application expected to be of the order of the 
typical hadronic mass scale (> 1 GeV), where the optimization is verified and, in addition, 
the prediction turns out to be reasonably stable. After the Borel transformation, possible 
subtraction terms are eliminated and eq.(EO) can be rewritten as 



fl(r,T',y) = J dvdv' e ^)p(u,u',y) . (3.9) 

Eq. ( |3.3| ) shows that the preliminary evaluation of the constants F(fi) and F + (fi) is 
necessary to exploit the sum rule for the determination of t\/i- This calculation is discussed 
in the next Section. 



IV. DETERMINATION OF F+(p) AND A+ 
The QCD sum rule determination of F + (fi) can be done by analyzing the correlator 

#(o/) =i fd 4 x e ik '- x < 0\T[J v s '{x)jf (0)^0 > , (4.1) 



whose dispersive representation takes contribution from the + qQ pole 

= 9A + ■ + - / M- + subtr. , 4.2 

2A+ — uj 1 — ie it Ju' c v — uj' — ie 

u' c being the effective threshold separating the contribution of the first resonance from the 
continuum. 
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It is straightforward to derive, at the next-to- leading order in a s , the contributions to 
the perturbative and non-perturbative parts of In the MS scheme one obtains 



-Im* per V) 

7T 



3oo 



■e(o0 



1 



2a s( //iv 17 , 2tt 2 



7T 
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and, considering non-perturbative vacuum condensates up to dimension five, 

#< GG >(u/) = 



7T 



^ — 

yr,<gGq>( /\ = m <qq> 
1 J " 2a;' 3 



(4.3) 

(4.4) 
(4.5) 
(4.6) 



where the relation < qg s a ■ Gq >= m 2 , < qq > has been used (mjj = 0.8 ± 0.2 GeV^ 2 0). 
Consistently with the first order in a s considered here, we neglect perturbative corrections 
to the coefficients of the higher-dimensional condensates in ( (4.5[) and ([4.6|) . The scale- 
dependence of the quark condensate is 



< qq > 



<">=<*> 



where d 



12 



33 - 2n/ 



a s (/x) 



v a s (/i ) 

is the number of "active" quarks and 
12- 1 L 153 -19n, ln ( ln (4))l 



(4.7) 



33 - 2n f l n [ (33 - 2n/) 2 l n 



(4.8) 



The numerical value we use for the quark condensate at fio = 1 GeV is < qq > (/io) 
(-240 Mel/) 3 . 

The Borel improved sum rule for F + reads 



[F+(/z)]V 



,A+ 



In 2 Jo 



di/ v 



I ..'2 



7T V K V' > 6 



17 



2tt' 



e r' 



+ < qq > (ti) 



i + 



2ct s 



rrir 



7T 



4r' 2 



(4.9) 



Eq.( [4.9|) shows that, for the renormalization scale /i, hence for the argument of a s , of the 
order of a typical strong interaction scale, ~ 1 — 2 GeV, the next-to-leading contribution 
to the perturbative part of the sum rule for F + is large, similar to the situation met in the 
case of F §-0. 

A /^-independent constant F+ n can be defined, using eq. ( |4.9| ) and the relation between 
F + (fi) and the matrix element of the scalar current in full QCD, in the same way one defines 



a \x- independent leptonic constant F ren for the sf = ~ doublet ]T2 
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Kin = [«.(aO]* [i - — z] (4.io) 

153 - 19n/ 381 - 30n/ + 28vr 2 



where, in the MS scheme, Z = 3 . 

(33 -2n/) 2 36(33 - 2n f ) 

From the above equations, a determination of F + (fi) and _F+ n can be obtained, together 
with the parameter A + . The latter quantity, for example, can be evaluated by considering 
the logarithmic derivative of eq. ( |4.9|) with respect to the Borel parameter r'. With r' 
in the range (1 — 2.5) GeV and the threshold v' c in the range (2 — 3) GeV, and choosing 
Aqcd = 380 MeV, we obtain for A + and F+ n the curves depicted in fig.|I|. The corresponding 
predictions are: 

A + = 1.0 ±0.1 GeV F+ n = 0.7 ± 0.2 GeV* . (4.11) 



The expressions relevant to sum rule for the constant F(/j) in (|3.4j) referring to the 
(0 - , l - ) heavy meson doublet, are given by eq. ( |4.3| ) for the perturbative part, and by 
reversing the signs of ( |4.4| - [4.6|) for the vacuum condensate contributions. An extensive 
analysis of this quantity can be found in || . We only repeat here the numerical calculation of 
0, using the same input parameters adopted for F + , and choosing the continuum threshold 
in the range v c = 2 — 3 GeV. The result is 

A = 0.5 ± 0.1 GeV F ren = 0.45 ± 0.05 GeV? . (4.12) 

The difference A = A + — A corresponds to the difference between and m£> , where D and 
Z) are the spin averaged states of the ~ and | + doublets. Our central value A = 0.5 GeV 
allows us to predict ~ 2.45 GeV with an uncertainty of about 0.15 GeV. 

It is worth reminding that determinations of the leptonic constant F + at the order a s = 
by QCD sum rules gave the result: F+ = 0.46 ±0.06 GeV§ f30fl and F + = 0.40 ±0.04 GeV§ 



and A+ = 1.05 ± 0.5 GeV or A+ = 0.90 ± 0.10 GeV H,|37| depending on the choice of 



the interpolating currents. The difference with respect to the values in ( [4.11| ) is the effect 



of sizeable radiative corrections. However, as far as the determination of Ti/ 2 is concerned, 
since radiative corrections affect both the three-point correlator, and the two-point functions 
determining the leptonic constants, it is still possible that a partial compensation occurs in 
the ratio determining the form factor; we shall see in the following that this is, indeed, the 
case. 

Determinations of F + by quark models [SB] give results in agreement with ([4=. 1 1|) when 



relativistic models [[39[] are employed: F + ~ 0.6 — 0.7 GeV 3 ^ 2 . Conversely, lower values are 
obtained: F + ~ 0.235 Gel/ 3 / 2 using non relativistic models [|40[j . 

The values in (|4 .11 ) and ( 4.12 ), or the equations corresponding to the respective sum 



rules, can be used as an input in (|3.3| ) to determine Ti/ 2 



V. 0(a s ) CORRECTIONS TO THE SUM RULE FOR r 1/2 

In order to calculate 0(a s ) corrections to the perturbative part of the sum rule for ri/ 2 (y), 
one has to compute the two-loop diagrams depicted in fig.@. Also the non perturbative term 
proportional to the quark condensate receives 0(a s ) corrections, as discussed in Section VI. 
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On the other hand, consistently with the order in g s considered here and with the previous 
estimates of F + and F, radiative corrections to the contributions of higher dimensional 
condensates will not be included. 

We start from the calculation of the perturbative part. As shown in [12], it is useful to 
directly deal with the double-Borel transformed expressions of the integrals corresponding 
to the various diagrams, a procedure which considerably simplifies the resulting calculations. 
This is the strategy we follow here, keeping different values of the Borel parameters r and 
t' in ( p.9| ). Adopting the standard dimensional regularization procedure, we compute the 
diagrams in D space-time dimensions, using the Feynman rules of HQET, and then we 
consider the expansion for e = (D — 4)/2 —>■ + . 

At the order a s = the expression for the Borel-transformed correlator (|3.1|) in the 
variables u,u' is given by 



D, 



iiv 



->- T -;, 

T 



AN C F(D/2) 



T 



D-l 



(4tt) d / 2 [V 2 {t/t')] d / 2 



(5.1) 



u 2 + 2uy + 1; r and r' are the Borel 



(|5.1|) shows that one has to perform the 



where N c is the number of colours and V 2 (u) 
variables related to uj and u', respectively. Eq. 
calculation with r ^ r' from the very beginning; a different situation is met in the case of 
the Isgur-Wise form factor £(?/), where the choice r = r' is allowed by the symmetry of the 
three-point correlator, with a remarkable simplification of the analysis. The requirement of 
keeping different values of the Borel variables represents the main technical difficulty in this 
calculation. 

In the following we give the results for the various diagrams in fig.|2], together with few 
details concerning the calculation; some useful formulae are collected in the Appendix. The 
overall computational strategy follows that adopted in for the calculation of the a s 
corrections to the Isgur-Wise function £(y), and we refer to this paper for further details. 



A. Diagrams D\, D2, D% 



Let us first consider the diagrams Di, D 2 and D 3 , where the gluon has both vertices on 
the heavy quark lines. Applying HQET Feynman rules we obtain for the the diagram Di. 



D l = mN c g 2 s C F yTr 



75" 



1+ 1 



-7^75- 



1+ 



-7 



(5.2) 



d u s d u t 



(2tt) d (2tt) d {uj + 2v ■ s)(u + 2v ■ t){u' + 2v> ■ s)(u' + 2v' ■ t)s 2 {s - t) 2 
and, after double-Borel transform in the variables u, uj', 

du /t\d-2 r 1 du ] 




D-A-' ^ 



V 2 (u^) 



D/2 



/T\D-2 r 



D/2-1 



(5.3) 



where 
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-2D-5 



z r = wiN^c^r^ - i)r(f) 



(4 



D/2 



(5.4) 



and CV : 
and D 3 : 



2N C 



. By a similar calculation, after borelization we obtain for the diagrams D 2 



D *-- Z ^ l - T -) (D-3KD-4) 



(5.5) 



3 

1=1 



„ / T\ 1 ,t'\D-4 

»\ t>) (£>-3)(£>-4) V T - 



(5.6) 



Expanding ([5.31)- (^75[) around e = 0, and summing up the contributions of Di — D 3 , we 
obtain: 



■(1-3)] -[vr{y)-i 



(5.7) 



^(3/' 7) + h {v> 77)] + ln ^ 2 (^:) r (2/) - y\fy 2 - lr2 (y) 



T 



-2yJy*-lF 1 - )F 1 - + 4 - In - + C?(e 



r\2 



where r(y) 



is the function already met in the expression of the Wilson coefficient 



in eq.(^.3|), and the variable y± is defined as y± = y ± \fy 2 — 1; the function /i(y, z) is 



%, 2;) 



L 2 (l-y 2 _) + L, 



y + (i-y 2 -) 
y+ + z 



The expression of the diagram D4, involving the light quark self-energy, is 



D 4 = UN c g 2 s C F {2 - D)Tr 
d D s d D t 



. i+ j i+f 
75^—7^75 2 7 7 P 7 7 



(2tt) d (2tt) d (w + 2v ■ s)(u' + 2v> ■ s){s 2 )H 2 {s - tf 



(5.8) 



with L 2 the dilogarithm function: L 2 {x) = — J X ^t 111 ^ 1 - . The function -Fi(z) in ( |5.7|) is 
reported in the Appendix. 

B. Diagram D4 



(5.9) 



Integration over t is straightforward. The double-Borel transformed D4 is immediately 
obtained: 
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t'> D-A 



which, after expanding around e = 0, becomes 



v2 {t) 



D/2-2 



(5.10) 



l_Il nF 2 (l)+0( e ) 
2e 2 V' v > 



(5.11) 



C. Diagrams D$,Dq 



The most difficult diagrams to compute are D 5 and D e , where the gluon has one vertex 
on the light quark line. Starting from the expression of -D5 



D 5 = 8iN c g 2 s C F Tr 



1+ 



1+ i' 



d D s d D t 



75^—7^75^- 



Sfyt 



(5.12) 



(2ir) D (2n) D (lu + 2v ■s)(w + 2v t)(u' + 2v' ■ s)s 2 t 2 (s - t) 2 
and using the identity 

1+ 1 



Tr 



■ 1+ i 1+ f a , 
75—^—7^75 2 7 VT 



Tr 



75 2 7/^75 2 



1 + // y /3 + -[ 7 a , 7 ^]+2A a 



(5.13) 



one can write: D5 = D5 + Z^ 2 " 1 + in correspondence to the above three terms. The 
first one, obtained by using s ■ t = |[s 2 + 1 2 — (s — t) 2 ], is given by 



(1) ^ n 1 



2 D-2 



r -) D - 2 V 2 (^ r 1 



T 



r> jo 



du 



T 



v^i) 

du 



1 v 2 ( T - 



D/2-1 D _ 3 > \ T 



T , \D-2 r _^ / -T\-\D/2 /•! (1 - m) 2 ^ 



D/2-1 



(5.14) 



where it is worth noticing the nontrivial dependence on the Borel parameters r and r'. Eq. 
( |5.14| ) can be simplified by integrating by parts the last integral 



i , (l-u) 2 - D 
du 



1 D - 2 t' 



o 



V\u^) 



D/2-1 



+ 



D- 3 D-3 



's /"I 



r v r / Jo 



y + —) du 



[l-u 



,3-D 



V 2 H) 



_D/2 



D-3V^ Jo 
and by considering that, for D = 4 + 2e, one can write 



D/2 



(5.15) 



13 



-i (1 - u) 4 ~ D 
du 



to 



V 2 (u^ 



D/2 



v 2 { — 



_ Y{Fl (L)_ e[2Fl (L) +F .(L)] } + 0(e 



■T 
T 



(5.16) 



with the functions F^iz), F&(z) and Fq(z) reported in the Appendix. Moreover, the first 
integral on the r.h.s. of eq. ( |5.15|) can be evaluated by performing one more integration by 
parts and making use of the identity 



,i ln w (l -u)(u^+y) 
D / du- 



r r 1 . In^Hl-u) 



D/2+1 



T JO 



n — / du 



l-u 



V 2 { T - 



D/2 



D/2 



(5.17) 



with the result 



i , (l-u) 3 - D 
du 



-1 



o 



V 2 [u T - 



D/2 



V' 



2 t 



2e 



F 9 -)+-lnV z - 



T 



T 



+ e 



W^(7) + <) lnVa (7) + ^(7) + ^(7 



+ 0(e) 



(5.18) 



(the expressions for the functions Fr(z),Fs(z) and Fg(z) can be found in the Appendix). 
The resulting expression for (the corresponding contribution to D G can obtained anal- 
ogously) appears rather simple, in spite of the involved expressions of the intermediate 
formulae; matter of fact, one has 



D 



(i) 



1 T 



T 



T 



T ' T 



0( t 



(5.19) 



,(2) 



The contribution D 5 , 

Df = 8tN c g 2 C F Tr 

d D s d D t 



75 



1+ "-- 1+ "'W] 



S a t/3 



2 2' 



(5.20) 



(2tt) d (2n) D (u + 2v ■ s)(cu + 2v ■ t)(u' + 2v' ■ s)s 2 t 2 (s - t) 2 ' 
after double-Borel transform, can be written is terms of a triple integral: 



D 



(2) 



HN c g 2 s C F (l + y)(v-v') 



rp-i) fVr Al /* 



(4tt 



dAi — / dz\ 
V Jo 



dz- 



Z\\Z X ZiZ 2 Z 2 



■ D/2-1 



(k 



2\D-1 



(5.21) 



where: Zi — 1 — z*; fc 2 = zi(g 2 Z2 +P 2 ^2) — ^i(?^2 +P-22) 2 , with p = XxV and g = - + p. The 
integral in (|5.21 ) can be simplified by changing the variables: 



Ai — - 

T 



Z\ 



^2-2(1-^ 

UiU 2 - 1 



Z2 



Ml 1 

A 1 



Au 2 -2(l-t) 



(5.22) 
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and by integrating over A, obtaining: 



D 



(2) 



Z» r (1 + y) 



V 2 - 



2 t' D - 3 

r(D - i) r 00 
r(f )r(f - i) A 

r\-D 1 



£>/2 



( Ul « 2 - l)^ 2 



-2 



D/2-1 



1-u 



D-l 



(5.23) 



D/2 



aw 



D/2 



The corresponding contribution to can be obtained analogously. Also in this case it is 
important to notice the dependence on the ratio of Borel parameters r and r'. 



(2) * (2) 

When we expand D\ and D\ in e we obtain a quite simple result, namely: 



D 



D 



(2) 



(2) 



r r 1 



2e 



T ' T 



T 



T 



0(e) 



1 +^ 2/r 



2e 2 vr 
and, for the sum ^=12(^5 + D^): 

E fcl , 2 (4" + cf) = f{(l-p)i 



0(< 



- (l"3)lnV 2 Q + (l + y)(l 



T 

T 7 



r v r 



Let us finally consider the third contribution to D§: 



£>( 3) = 8iN c g 2 s C F Tr 



1+ ^ 1+ ^' 
75^—7^75 — 



-r 



t' V 



"In 



0( t 



where 



s„ 2t> • t 



(2vr) D (2tt) d (w + 2v ■ s)(uj + 2v ■ t)(u' + 2v' ■ s) s 2 t 2 (s - t) 2 
The integral (gjg) can be related to simpler ones by using the identity 



a, /" d s a t 



+ 



(5.24) 
(5.25) 



(5.26) 



(5.27) 



(5.28) 



(2tt) d (2tv) d (lu + 2v s)(oj' + 2v' ■ s) s 2 t 2 (s - t) 2 

d D s d D t 2vts a 
(2%) D (2tt) d (lu + 2v ■ s)(u + 2v ■ t){u' + 2v' ■ s) t 4 (s - tf 

d D s d D t 2v ■ t s a 

(2n) D (2vr) D (u + 2v ■ s)(u + 2v ■ t){u' + 2v' ■ s) s 2 t A 
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d D s d D t 



(2tt) d (2ir) D (lu + 2v t) 2 (u' + 2v' ■ s) s 2 t 2 (s - t) 2 

= 4 1) +4 2) + 4 3) + ^ ) - 

After Borel transformation, the results for jW, i — 1, 2, 3 are given by 
4 t 2D - 5 /D x /Z)w r a 1 

r(--iW_Vv + -y 



(D-2)(D-4) (4tt) d V2 



7-' /a 



D-2 



_2D-5 



D 



V 2 r 



D/2-1 ( 4?r )D V 2 

+ V) 



r(-- i )J (- 



D 



(5.29) 



(5.30) 



(5.31) 



7-' / Q 



D-2 



t'\D-2 f 1 



du 



D/2-1 



a 



V2 f T _,yw* + ( 1+ y^)h) ' ~<L du rw t>\~\ d / 2 



V 2 (u^) 



Ji 3) = 2 



T 



2D-5 



(4tt) c 
_/ 



^t-'Mt) (7) "«/„ du 



3-D 



D/2 



(5.32) 



3-D 



V 2 {u T - 



D/2 



As for JW, it can be obtained as the result of a differential equation introduced in f4"2"| , |T2| 
which gives 



D-2 (4tt) d V 2 

9 s " 



2 

D-2 



(5.33) 



v2 (t) 



D/2-1 D _ 4 



TX4-D 



(v + ^t/ 



-(D-4) 



— ) / d'i 



u 



D-4 



D-2 /r\2D-5 r 1 



D/2-1 



D-4 



u 



D-4 



D/2 



r' \ 4-D 



It- 



(«' + 



r 
u — v 

f / a 



Summing up the four contributions one can write 
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D 



(3) 



D-4 



where 



D 



(3.1) 



1 



(D-2)(D-4) 



D/2-2 



T 



(5.34) 



(5.35) 



D 



(3.2) 



^ - 



D-2 



t'\D-2 r 1 



+ 



(l + 2/-)(-J / duy^- ~~i 



du- 



D/2-1 



D 



(3.3) 



D 



(3.4) 



T \D-1 
T 

1 







D/2 /•! (1 — It 

a« 



.3-D 



r \3-D 
D-2 I V 



r 



r 



D-2 



+ 



D-4 



D/2 V 1 " U 7 



r 



1 



(5.36) 



(5.37) 



(5.38) 



-(D-4)( 







D/2 



















dtt- 



D-4 



D/2-1 



D-2 r 1 , it 

aw — 



D-4 Jo 



D-4 



l/ 2 U r - 



D/2 



T\4-D 
U — — 1 

T 1 



1 — 11- 



T 



The corresponding contribution to Dq is obtained analogously. It is worth noticing that 
some of the integrals in D® and Dg 3 "* require an expansion up to order e 2 to take care of 
the D — 4 factor appearing on the l.h.s. of ( |5.29 ). 

~ (3) " (3^ 

Despite the involved structure of the expressions for the various terms in D\ and D\ , 
the result for the sum Dg + D 6 is rather simple: 



!-3 +(1+2/) 1 + 3 -7^3 



r\rr „ / r\ r' /r' 
— iM — 

T' J L T V T' ' T v T 



+ (i - T ~) [% 2 - ^(^(7) + 1- 2 + 2 + 21n ^)] + . ( 5 - 39 ) 



1 



where one can notice that spurious — terms cancel out 



D. Final result 

The sum of the contributions of the diagrams D\ — D 6 gives the result: 
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*> = t A = ^ (i - £){ - \[\ - vr(y)} + (i + »)(i + $)g(L 

i=l \ 



+ 



h {y^) + h (y^)} + y r (y) ln ^ 2 ( 7) - v\lv 2 - lr2 (y) - ln (7) 2 



-7T 



+ 0(e) 



(5.40) 

where = (xF%(x) — ^Fi(-))/ (1 — x). The important point to notice is the structure of 
the 1/e singularity in ( p.40[) , which does not depend on the Borel parameters t,t'. 

Eq. (|5.40| ) represents, for all values of the Borel parameters r and r' the 0(a s ) correction 
to the triangle diagram representing the correlator ( |3.1| ). In particular, in the limit r ~ r' 
the expression can be simplified, giving 



a* 



D = D — - 

7T 



H 


n i(y)- 


+ 2 7£ ; 


N 7(2/)- 


+ 2.„(^) 


h 7(2/)] 




2 - 


2 - 




L 2 J 



+ ^yh{y) 



+ 



|(y 2 - 1) - yJy^T\r 2 (y) + (|yr(y) + ±) In[2(l + y)] + ^vr 2 + || (5.41) 

where 7^ is the Euler constant and i = - + 7^ — ln47r; 7(1/) was defined in (|2.5|) , and 
/i(y) = h(y, 1). In theMS 1 subtraction scheme the 4 pole cancels with the renormalization 
factors of the heavy- light and heavy-heavy quark currents in the correlator ( |3.2|) [p,|43",44|: 
Z\ x = 1 — % , iT^ = 1 + f^j(y), so that the finite part represents the correction to 
the Borel-transformed correlator we are looking for. 

It is possible from eq.( |5.41| ) to determine the spectral function fp ert in ( |3.7| ) at the order 
a s , which is required to perform the continuum subtraction in the QCD sum rule analysis. 



After changing the variables in (|3.9| ) to a± — v ± z/, and integrating in er_, with integration 
limits < cr + < +00 and —ra + < cr_ < ra + 

function proportional to a\ pi(y) + P2(y) m \-^~ 



y+l 

so that 



^tt), one is left, for r ~ r', with a 



D oc pi(y) + 



11 

IE 





(5.42) 



Comparing (|5.42|) with (|5.41 ) we get, with Dq the a s = term flBTI 



Pi (2/) 



a* 14 



A)^ 3^(2/) + [g(r - 1) - yyjy 2 - l\r\y) + (-yr(y) + -) ln[2(l + j/)] 



8 11 



7(2/) 



p 2 (y) = 2D t 



a, r 



7T 



1 - 



3 

7(2/) 



(5.43) 
(5.44) 



and therefore the spectral function p veTt {y, v\ y) 
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VI. NON-PERTURBATIVE CONTRIBUTIONS AND FINAL SUM RULE 



Once the perturbative contribution to the sum rule has been computed, one has to 
consider non-perturbative power corrections, and, as anticipated, we include the vacuum 
condensates up to dimension five. 

The lowest dimensional term in the expansion of U np in eq. ( |3.7| ) is the quark condensate 
< qq >. At the tree level, the Borel transformed result for this contribution is simply given 
by 



n<gq> 



-i < qq > (v — v 



i' ■ 



(6.1) 



The 0(a s ) correction to ( |6.1| ) is computed from eight diagrams obtained from those in 
fig. 2 replacing the light quark line by the quark condensate contribution to the relevant 
propagator. The correction reads 



J2 D 

i=l 



<qq> 



Dr q> {-)H{ry) 

v 7T ' 



(6.2) 



where 



H(t, t') = - - 2[1 - yr{y)] lE - [1 - yr(y)} (in (£) + (£)) + In (1) [1 + yr(y)} 



+ \nV 2 (^) -5 + (l + y)r(y)+y -h(y,^) - h(y, ^) - r(y) ^V 2 {^-) 



T 
T 



'T 
T 



(6.3) 



with the notations previously defined. The calculation of the dimension five contribution is 
straightforward. Then, the final sum rule for n/2 is 



2A 2A ' / 

2 r 1/2 (y,/i)F(^F + (//) e""-— = ^ dvdv' (P^\v, v\ y)e 



v v 

T T > 



<qq> (/i)(l + — H(r,r 



j\ mg r 1 1 4y + l . 
2 l 2r 2 2r' 2 3rr' 



(6.4) 



where 



v — V 



16tt 2 (y- l)v^^T 



1 + - ((2 - 7 (y)) In + ^ + ^ + c p (y))] 8(1/ - i/_ ) 0(i/ + - i/) 



(6.5) 



and 
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Cp(y) = ^yh(y) + (Jo 2 - 1) - y^y 2 - l)r 2 (y) 



+ ~ y + ln(2(l + y))) + ln(2(l + y)) (6.6) 

(f± = y±v)- The integration domain .D is constrained by the conditions v < u c , v' < u' c . 

Since the form factor Tin is defined by the matrix elements of weak currents in the effec- 
tive theory, it depends on the subtraction scale /i, and the sum rule (|6.4|) clearly reproduces 
this feature. As discussed in Sect. II, and in analogy with the case of the Isgur-Wise func- 
tion £, it is possible to remove the scale-dependence by compensating it by the analogous 
//-dependence of the Wilson coefficients relating the b — > c axial current in full QCD to the 
dimension 3 currents in HQET and by defining rjyJJ as i n eq.( |2.7| ). This is the function we 
shall consider in our numerical analysis. 

VII. NUMERICAL RESULTS 

The numerical analysis of the sum rule for t\i 2 can be carried out using the same input 
parameters adopted in Sect. IV for the determination of F + . In particular, we use the explicit 
expressions of the two-point sum rules determining the leptonic constants F + and F that 
appear in the pole contribution of eq. (|3.3| ). We vary the threshold parameters in the ranges 
u c = 2 — 3 GeV and v' c = 2.5 — 3.5 GeV, obtaining an acceptable stability window, where the 
results do not appreciably depend on the Borel parameters, in the ranges around r ~ 1.5 
GeV and r' ~ 2 GeV, respectively. The contribution of the nonperturbative term in the 
three-point correlator represents a small fraction of the total contribution; on the other 
hand, the a s correction in the perturbative term is sizeable, but it turns out to be partially 
compensated by the analogous correction in the leptonic constants F and F + . Notice that 
this is a remarkable result, not expected a priori since the normalization of the form factor, 
for example at zero recoil, is not fixed by symmetry arguments. The perturbative corrections, 
however, do not equally affect the form factor for all values of the variable y, but they are 
sensibly y— dependent, with the effect of increasing the slope of Ty 2 with respect to the case 
where they are omitted. 

The results for r^iy) are shown in fig. 3, where the curves refer to various choices 
for the continuum thresholds. The region limited by the curves essentially determines the 
theoretical accuracy allowed by the present calculation. 

Considering the y dependence, the limited range of values allowed by the mass difference 
between D and Dq permits the expansion near y = 1: 

r^(y)=T 1/2 (l)(l-pl /2 (y-l)+c 1/2 (y-l) 2 ) . (7.1) 

A two-parameter fit to fig. 3, in terms of the normalization at zero recoil and the slope, gives 
r i/2(l) = 0.31 ± 0.06 and p\j 2 = 1.5 ± 0.4. The inclusion of the quadratic term modifies the 
fit as follows: 

r 1/2 (l) = 0.35 ±0.08 , p\ /2 = 2.5 ± 1.0 , c 1/2 = 3 ± 3 (7.2) 
which is the result we quote for our analysis. 
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The immediate application of this result concerns the prediction of the semileptonic B 
decay rates to D and D\. Using V c b = 3.9 x 10~ 2 and t(B) = 1.56 x 1CU 12 sec, we obtain 

B(B -> D £u) = (5 ± 3) x 1(T 4 B(B -> D\lv) = (7 ± 5) x 10~ 4 . (7.3) 

This means that only a very small fraction of the semileptonic B — > X c decays is represented 
by transitions into the sf = ^ + charmed doublet. Although small, however, one cannot 
exclude that such processes will be identified, mainly at dedicated 5-facilities which will be 
running in the near future. At present, the measurements of semileptonic B — ► D** decays 
only provide data on the members of the sf = | + doublet [fl5|Jl6[j , since the doublet with 
sf = i + is not distinguished from the non-resonant charmed background. In particular, in 
f45|| the B semileptonic branching fraction to the final states Dtt and D*n of (20 ± 5) 10~ 2 
is reported. 



VIII. CONCLUSIONS 



We conclude observing that HQET has proven to be a powerful tool to handle heavy 
quark physics. However, predictions derived in this framework should always be supported 
by the computation of 1/ vtlq as well as radiative corrections. The role of both depend on 
the specific situation one is facing with. For example, they turn out to be important for the 
B meson leptonic constant, while they are moderate for the Isgur-Wise function, as derived 
in ||12j| . We have presented here the case of the universal form factor Tyz{y) describing 
B semileptonic transitions to the excited J p = (0 + , 1 + ) charmed states, using QCD sum 
rules in the framework of HQET. As already shown in ||12|| , the computation of loop integrals 
results to be greatly simplified within HQET. The task of computing perturbative corrections 
to Ti/ 2 (y) is justified by manifold interesting phenomenological features of orbitally excited 
states as well as by the many theoretical interests already mentioned. We have obtained a 
situation similar to the case of the Isgur-Wise function, namely radiative corrections are quite 
under control for 71/2(2/), while they affect considerably the value of the leptonic constant 
F + of the s t = 1/2 doublet. 
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APPENDIX A: PARAMETRIC INTEGRALS 



The calculation of the two- loop diagrams relevant for the form factor t\ii essentially 
follows the analogous calculation for the Isgur-Wise function [12], with the main difference 
represented by the need of keeping different Borel parameters, due to the non-symmetric 
nature of the problem at hand. We only recall here that, in momentum space, the Feynman 
rules of HQET can be reduced to the heavy quark propagator ^rr-x, and to the 

heavy-quark-gluon vertex ig s v^; N c is the number of colours and TV[4f 4r] = Cp = J 1'C 1 . 

The calculation of the loop integrals is performed in D = 4 + 2e Euclidean space-time 
dimensions. The main ingredients are the representations of the propagators of the massless 
quark and of the heavy quark: 
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(ve, se obtained from the four- vectors v , s by a Wick rotation); in particular, 
for the computation of the integrals after Borel transformation, since 
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(A3) 



The master integrals needed in the evaluation of the loop integrals can be found in |f47U44| . 
Here we report a number of parametric integrals useful for the calculation of the Borel- 
transformed expressions D\ — Dq: 
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In Fi — F 8 the combinations have been introduced: 
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where L 2 (x) is the dilogarithm. The integrals i*i(/z) — Fg(z) coincide with those reported in 
12]| for z = 1. Finally, a useful identity needed in the calculation of , D 6 is 
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FIGURE CAPTIONS 



Fig. 1 

Binding energy parameter A + and leptonic constant F+ n of the doublet sf — \ , from 
the QCD sum rule analysis of the correlator eq.( |4.1| ). The curves refer to three choices of 
the threshold parameter v' c : v' c = 2.5 GeV (continuous line), v' c = 3.0 GeV (dashed line), 
v' c = 3.5 GeV (dotted line). 

Fig. 2 

Two-loop diagrams relevant for the calculation of 0(a s ) corrections to the perturbative part 
of the QCD sum rule for the form factor ri/2- The heavy lines represent the heavy quark 
propagators in HQET. 

Fig. 3 

The universal form factor T™%(y). The curves refer to choices of the threshold parameters: 
v c = 2.0 GeV, v' c = 2.5 GeV (continuous line), v c = 2.5 GeV, v' c = 3.0 GeV (dashed line), 
v c = 3.0 GeV, u' c = 3.5 GeV (dotted line). 
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